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By exact diagonalization, we investigate several spin-i 
Ji — J2 real clusters of molecular scale with different shapes. 
Our calculations show that when the ratio, 77, of next nearest 
neighbor to nearest neighbor bonds is equal to 1, even the 
cluster of only 25 sites exhibits the bulk behaviors and has 
the quantum phase transitions. Two effective critical points 
are around J2/J1 = 0.3762 and 0.612 respectively. They 
are very close to those of the infinite Ji — J2 square lattice. 
But, when r\ < 0.85, the quantum phase transition around 
J2/J1 = 0.3762 disappears. By calculating the distributions 
of the average values of Sf , the bulk behaviors are demon- 
strated graphically. In the intermediate phase, the sites on 
the corners have distinctly different character from the other 
sites. The distribution is obviously centralized on the corner 
sites. 



The two-dimensional spin-i j l j 2 Heisenberg model 
has been the object of intense investigation through years 
since the antiferromagnetism can be destabilized by frus- 
trations and it exhibits rich physics and diverse proper- 
ties. Its Hamiltonian reads 



H 



Ji Si ■ Sj + J2 ■ Sj> 



(i) 



(n.n.) 



{n.n.n.) 



where (n.n.) and (n.n.n.) are the sum over the nearest 
neighbors and next nearest neighbors respectively, J\ > 
and J 2 > 0. For the two-dimensional antiferromagnet 
without frustration (J 2 = 0), notwithstanding a mathe- 
matically rigorous solution is lacking, it is by now well 
established that the ground state (GS) has Neel order 
by many works |]|l5|]. When J2 > 0, the competition of 
purely magnetic interactions can lead to the destruction 
of the long-range order. Many analytical approximations 
|p]-^,P|-[T3| and numerical techniques (|-||] are applied to 
investigate the phase diagram of the J\ — J2 model. For 
small J2/J1, the GS will keep Neel order described by a 
wave vector (tt, tt) since the frustrations are too weak to 
destroy it. When J2/J1 is large enough, the GS is dom- 
inated by the n.n.n. interactions and has a certain type 
of collinear order described by (tt,0) or (0,7r). But, for 
intermediate J2 j J\ , the competition between Ji and J2 
terms results in a spin-liquid GS with a gap to magnetic 
excitations. The two critical points of phase transition 
are at J2/J1 = 0.38 and 0.60 respectively 



The situation of real cluster, when it has small size, 
may be different from that of infinite system (the word 
'real' means that the free boundary conditions are taken). 
Especially for real clusters of molecular scale , which has 
only tens of sites, the boundary effect might be heavy 
since many sites are on its boundaries. For instance, the 
4 x 4, 5 x 5 and 6x6 real clusters have 75%, 64% and 
56% sites on their boundaries respectively. In order to 
understand the spin status of real clusters, we introduce 



<3(?1 



Si ■ Si 



G)exp(^-(r--f i )), (2) 



as a measure of magnetic orders, where TV is the total 
number of sites, fi and fj the coordinates of sites i and j 
respectively, \G) represents the GS. At first, we calculate 
Q(ir,ir) of the 5x5 real cluster by exact diagonaliza- 
tion. Our calculation shows that Q(tt, tt) decays continu- 
ously and monotonically in the region of J2/ J\ ~ — 1.0 
(curve (e) in Fig. 2). Especially, unlike the case of the 
infinite system, the real cluster has no quantum phase 
transition around J%j ' J\ = 0.38. One may think that 
the cause is too small freedom to form the bulk behav- 
ior of spins. It is meaningful to investigate whether the 
real clusters of molecular scale always keep this property. 
However, the real cluster with other shape might have 
different properties due to the boundary effects. If de- 
fine n = N* n n n d /N* n nd , where N^ nd and N*° n n r ? represent 
the numbers of n.n. and n.n.n. bonds, we notice that 
r/ = 1 for the infinite J± — J2 square lattice, n may be 
an important factor for spin status. In other words, if 
a real cluster posses a structure of 77 = 1, it has possi- 
bly the similar quantum transition as the infinite system. 
When its boundaries are set up along diagonals, the real 
cluster is able to hold 77 = 1. Calculating Q(ir,ir) of 
cluster (a) (Fig. 1), we find that it is quite different 
from the real 5x5 cluster and Q(ir,ir) has a step-like 
change around J2/J1 — 0.3762 (Fig. 2). In a very small 
interval of J2/J1, which is less than 4 x 10~ 4 , Q(tt,tt) 
falls 55%. The abrupt change in such a narrow parame- 
ter region means that the small real cluster has possibly 
the bulk behaviors, i.e., there exists a quantum phase 
transition around J2/J1 = 0.3762. Larger Q(tt,tt) for 
J2/J1 < 0.3760 means the Neel-like order Jlq ], whereas 
small Q(ir,ir) for J2/J1 > 0.3764 means its losing. We 
will give the further evidence for this conclusion. Un- 



1 



like the infinite system has Q(ir, it) = when the AF 
long-range order vanishes, here zero Q(tt, tt) can not be 
achieved since the real cluster is finite. In the region 
0.3764 < J2/J1 < 1.0, although Q(tt, n) decays mono- 
tonically as J2/J1 increasing, it keeps positive. The real 
frustrated cluster always has a little antiferromagnetism. 

Generally, the bulk behavior occurs for these clusters 
with size large enough. The too small size may lead 
to its disappearance due to the limitation of freedom. 
But, our calculations show that when the real cluster 
of molecular scale keeps 77 = 1, it has the bulk prop- 
erties, especially the quantum phase transition around 
J2/J1 — 0.38. Even for the real cluster of 13 sites with 
the similar shape, there also exists a step-like change of 
Q(ir,ir) within 0.3381 < J 2 /Ji < 0.3886. It is noticed 
that cluster (e), i.e. the 5x5 real cluster, has 77 = 0.80. 
Three other real clusters (b), (c) and (d) (Fig. 1), which 
have the same number of sites but the different values 
of 77 (77 = 0.9, 0.89 and 0.85 respectively), are investi- 
gated. As 77 decreasing, the jump around J2/J1 = 0.38 
becomes smaller and smaller (Fig. 2). When 77 < 0.85, 
the step-like change disappears. Our calculations sup- 
port the conclusion that 77 is an important factor for the 
quantum phase transitions of real clusters and the en- 
hance of 77 advantages the bulk behavior. 

For the frustrated Heisenberg antiferromagnets, when 
the effect of J2 terms exceeds that of J\ terms, the n.n.n. 
interactions may result in the establishing of magnetic or- 
ders in its subsystems. If carefully checking the Q(w,w) 
curves of these five clusters, one can find inflexions at 
larger J2/J1 (Fig. 2). The difference of their geometric 
symmetries causes the critical J2 / J\ values obviously dif- 
ferent since for such small clusters, the symmetries may 
affect heavily the spin-Peierls states of the disordered 
phase. One notices that same as the infinite square lat- 
tice, cluster (a) keeps the rotation symmetry and 7/ = 1. 
Its singular point is at J2/J1 = 0.612(±0.002). The slope 
of Q(tt,tt) increases when J2/J1 < 0.612 while deceases 
when J2/J1 > 0.612. For convenience, the real cluster 
is divided into two subsystems A and B (Fig. 1). To 
understand whether there is quantum transition around 
J2/J1 w 0.612, it is useful to calculate Q(tt, 0) by formula 
(||). As shown in Fig. 2, Q(ir, 0) changes rapidly round 
J2/J1 ~ 0.612, and it keeps a obviously larger value for 
large J2/J1 than that for small J2/J1. It implies the ex- 
istence of the Neel-like order in both of subsystems A 
and B while no Neel-like order for the whole real cluster. 
When J2/J1 < 0.612, the competition of J\ and J 2 leads 
to the disordered GS. Namely, the GS has the character 
of spin liquid. 

The total spin of the GS of cluster (a) is nonzero since it 
has odd sites. It provide us another way to investigate the 
spin status. The spin-spin correlations lead to the proba- 
bility of spin up unequal to that of spin down on each site 
in the M space, here M is the magnetic quantum num- 
ber. Consequently Sf = (G + \Sf\G + ) ^ 0, here |G+) 



involves only the ground states in the M > subspaces. 
(It is quite different from the case of the GS with zero 
total spin, in which the spin up-down symmetry makes 
Sf being zero exactly on each site.) The distribution of 
Sf on sites depends on the spin-spin correlations. The 
study of it is helpful to understand the spin status. Fig. 
3 shows the cases of cluster (a) for J2/J1 = 0.37, 0.40, 
0.60 and 0.88 respectively. At J2/J1 = 0.37, the average 
value of Sf for the sites in subsystem A is 2.04 x 10 _1 
and its statistical fluctuation 5.48 x 10~ 3 , and that for 
the sites in subsystem B is —1.41 x 10 _1 and its fluctu- 
ation 1.43 x 10 -3 . i.e., the distribution of Sf is approx- 
imately uniform. It gives a good picture of the Neel-like 
order. But, when J2/J1 > 0.3764, the distribution is 
rather different. Fig. 3(b) and 3(c) shows the cases of 
J2/J1 = 0.40 and 0.60 respectively. For convenience to 
describe, we call the sites on four corners as corner sites, 
and the other 21 sites as the inner sites (Fig. 1). At 
J2/J1 = 0.40 and 0.60, Sf on each corner site is equal to 
1.88 x 10 _1 and 1.00 x 10 _1 respectively, whereas Sf on 
inner sites are distinctly small and equal to 2.63 x 10~ 2 
and 2.94 x 10~ 2 respectively. The fact that the magni- 
tudes of Sf on inner sites for J2/J1 > 0.3764 are much 
smaller than those for J 2 /.Ji < 0.3760 demonstrates the 
disappearance of the Neel-like order. For intermediate 
J2/J1, on each corner site, Sf always keeps positive, and 
its magnitude is much larger than those of the inner sites. 
Such distribution results from spin disorder (spin liquid) 
behavior. In this case, the spins on corner sites are in a 
special status. 

If J 2 / ' J\ increasing further, the n.n.n. interactions lead 
to the occurrence of a new ordered quantum phase. Fig. 
3(d) shows the case of J2I J\ — 0.88. In subsystem B, the 
average value of positive Sf is 1.93 x 10 _1 , and the fluc- 
tuation 4.67 x 10~ 3 . The four corner sites with negative 
Sf have the same value —1.37 x 10 _1 due to the geom- 
etry of cluster (a). The distribution of Sf in subsystem 
B exhibits the AF Neel-like order. In subsystem A, Sf 
are very small and their average value is 5.11 x 10~ 3 . It 
does not means that the spin status is disordered. By 
contrary, it is also ordered by the calculation of Q(ir, 0). 
A simple explanation to the smallness of Sf in subsystem 
A is given: cluster (a) consists of two subsystems A and 
B by Ji coupling. It is noticed that subsystem B has 
odd sites, but subsystem A has even sites. If J\ = 0, the 
distribution of Sf in subsystem B can be used to show 
the Neel picture from J 2 interactions. But, in subsys- 
tem A, we can not do it in the same way since Sf = 
exactly. For large J2/J1, the n.n.n. interactions domi- 
nate the physics of the cluster. Consequently, the site in 
subsystem A has small Sf. 

From our calculations, the following conclusions can 
be obtained. 1). although it has only 25 sites and nearly 
half of all sites on its boundaries, the molecular cluster 
(a) exhibits the obvious bulk behaviors. For small and 
large J2/J1, cluster (a) is in two different magnetic or- 
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dered phases. For the intermediate J2/J1, the spin status 
is disordered. Two kinds of quantum phase transitions 
take place around J2/J1 = 0.3762 and 0.612. respec- 
tively. 2). for cluster (a) with intermediate J2/J1, the 
spins on corner sites are in a special state. And their 
magnitudes of S| are remarkably larger than those on 
inner sites. We think that the larger real cluster with 
the similar shape as cluster (a) will keep this charac- 
ter. 3). 77 plays an important role for real cluster to 
keep the bulk behaviors. Its increasing advantages the 
formation of bulk behaviors. From our calculations, we 
think that only for real clusters with 77 near 1, it is pos- 
sible that there exists quantum phase transition around 
J2/J1 = 0.38. 

It is somewhat surprise that the two critical values 
J2/J1 = 0.3762 and 0.612 of cluster (a) are very close 
to those of infinite square lattice J2/J1 = 0.38 and 0.60 
|jl4j| . The fact that the infinite J\ — J2 square lattice 
has the same 77 = 1 as cluster (a) may be responsible for 
the result. 77 is a key factor for real cluster to possess 
the similar phase diagram as that in infinite system. By 
using exact diagonalization, H. J. Schulz et al. studied 
the 6x6 cluster under the periodic boundary conditions. 
Their calculations show that there is no step-like change 
of Q(7r,7r) around J2/J1 = 0.38, and Q(ir,ir) becomes 
obviously small only when J2IJ1 > 0.6 (see Fig. 3 (a) 
in Ref. ||). To obtain reliable critical values of infinite 
system, one must resort to scaling analyses like doing 
in Ref. ||. In addition, the 6x6 cluster has 56% sites 
on its boundaries. Although cluster (a) has remarkably 
fewer sites, the boundary effect might be weaker than the 
6x6 cluster since it only has 46% sites on its boundaries. 
The real clusters with such shape as cluster (a) might be 
better choice than the conventional n x n cluster under 
periodic conditions to study some properties of the infi- 
nite Ji — J2 square lattice. But, the calculable sizes of the 
clusters with such geometric shape are only 13 and 25. It 
is difficult to obtain the results of infinite system by ordi- 
nary extrapolation since the scarcity of available points. 
One notices that the region (0.3760 ~ 0.3764) of J2/J1 
of the real cluster of 25 sites, in which a step-like jump of 
Q(ir, 7r) takes place, is within that of the real cluster of 13 
sites (0.3381 ~ 0.3886). We think that the J 2 /Ji region 
of the larger real cluster will fall into 0.3760 ~ 0.3764. 
Then, we speculate that the critical value of the infinite 
system is J2/J1 = 0.3762(±0.0002). It may be more ac- 
curate than the previous result J2/J1 = 0.38. 

The part of calculations in this work have been done 
on the SGI Origin 2000 in the Group of Computational 
Condensed Matter Physics, National Laboratory of Solid 
State Microstructures, Nanjing University. This work 
was partially supported by the Ministry of Science and 
Technology of China under Grant No. |nkbrsf-gl999064|6. 
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FIG. 1. Five kinds of real clusters of 25 spins with different 
geometric shapes. All the sites except (a) 1 , 2 , 3 and 4 , 
(b) 1, 2', 3' and 4 , (c) 1, 2, 3' and 4, (d) 1, 2, 3 and 4 , 
(e) 1, 2, 3 and 4 have 1/2 spin. The solid lines represent the 
n.n. interactions Ji, and the dot lines the n.n.n. interactions 
J2. The cluster can be divided into two subsystems A and B. 
The solid circles denote the sites in A, and the solid and open 
squares those in B. 
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FIG. 2. Q(tt,tt) and Q(tt,0) vs. J 2 /Jl Curves (a), (b), 
(c), (d) and (e) are <3(7r,7r) for real clusters (a), (b), (c), (d) 
and (e) respectively (Fig. 1). Curve Q(tt, 0) is for real cluster 
(a). 



FIG. 3. Several kinds of Sf distributions on sites of cluster 
(a) (Fig. 1). The lengthes of arrows represent the magnitudes 
of Sf. a)Sf = 0^236 for the corner sites, b) ~Sf = 0.188 for the 
corner sites, c) Sf = 0.100 for the corner sites, d) Sf = 0.174 
for the center site. 
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